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The stability and collapse of fundamental unstaggered bright solitons in the discrete Schrodinger 
equation with the nonpolynomial on-site nonlinearity, which models a nearly one-dimensional Bose- 
Einstein condensate trapped in a deep optical lattice, are studied in the presence of the long- 
range dipole-dipole (DD) interactions. The cases of both attractive and repulsive contact and DD 
interaction are considered. The results are summarized in the form of stability /collapse diagrams 
in the parametric space of the model, which demonstrate that the the attractive DD interactions 
stabilize the solitons and help to prevent the collapse. Mobility of the discrete solitons is briefly 
considered too. 
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I. INTRODUCTION 

It is well established that the mean-field description 
of Bose-Einstein condensates (BECs) trapped in a deep 
optical lattice can be reduced, starting from the one- 
dimensional (ID) Gross-Pitaevskii equation with the cu- 
bic nonlinearity, to the discrete nonlinear Schrodinger 
(DNLS) equation [U, However, while the DNLS 

equation with the cubic on-site nonlinearity readily pre- 
dicts solitons it cannot describe the dynamical col- 
lapse, which is often observed experimentally in the self- 
attractive BEG H, 0. 

In the general case, the reduction of the 3D Gross- 
Pitaevskii equation for the BEG trapped in a "cigar- 
shaped" configuration to the ID form leads, instead of 
the "naive" cubic Schrodinger equation, to ones with a 
nonpolynomial nonlinearity p,|8|. Such models are useful 
in various settings, making it possible to obtain results 
in a relatively simple form, which compare well to direct 
simulations of the underlying 3D equations 01 ■ In partic- 
ular, the version of this model corresponding to the com- 
bination of the tight transverse trap and strong optical- 
lattice potential acting in the longitudinal direction takes 
the form of the discrete nonpolynomial Schrodinger equa- 
tion (DNPSE), which was introduced recently An 
essential property of both the continual nonpolynomial 
equation and its discrete counterpart is that they make 
it possible to model the collapse. 

A new variety of the BEG dynamics, dominated by 
long-range (nonlocal) interactions, occurs in dipolar con- 
densates, which can be composed of magnetically polar- 
ized ^^Gr atoms, as demonstrated in a series of experi- 
mental works |li|. In particular, the dipole-dipole (DD) 
attraction in the condensate may give rise to a specific 
d-wave mode of the collapse [l2|. On the other hand, the 
^^Gr condensate can be efficiently stabilized against the 
collapse, adjusting the scattering length of the contact in- 
teraction by means of the Feshbach-resonance (FR) tech- 
nique fllj. Another theoretically analyzed possibility is 



to create a condensate dominated by DD interactions 
between electric dipole moments induced in atoms by a 
strong external dc electric field [l^. A similar situation 
may be realized in BEG formed by dipolar molecules • 
In particular, recent experimental work TE\ has reported 
the creation of LiGs dipolar molecules in a mixed ultra- 
cold gas. 

A possibility of making 2D solitons in dipolar con- 
densates was predicted in several works. In particular, 
isotropic solitons and vortices [l^ may exist if the 
sign of the DD interaction is inverted by means of rapid 
rotation of the dipoles [l^. On the other hand, stable 
anisotropic solitons can be supported by the ordinary DD 
interaction, if the dipoles are polarized in the 2D plane 
[Tot . Solitons supported by nonlocal interactions were 
also predicted and realized in optics, making use of the 
thermal nonlinearity [2(3| . 

A natural extension of the consideration of the dipo- 
lar BEG includes a strong optical lattice potential, which 
leads to the discrete model with the long-range DD in- 
teractions between lattice sites [2l|, [g^]. In particular, 
ID unstaggered solitons in this model with the cubic 
on-site nonlinearity were studied in Ref. (2l| . It was 
shown that the DD interactions might enhance the soli- 
tons' stability. This conclusion suggests a possibility of 
suppressing the collapse by means of the long-range DD 
forces, but the study of the collapse is not possible in the 
discrete system with the cubic nonlinearity. The objec- 
tive of the present work is to introduce the DD interac- 
tion into DNPSE model and analyze the effects of the 
long-range interactions between dipoles localized at site 
of the discrete lattice on the soliton's collapse in the suf- 
ficiently dense self-attractive condensate trapped in the 
deep optical lattice potential. We here focus on unstag- 
gered solitons, leaving the consideration of staggered ones 
for another work. 

The paper is structured as follows. The discrete model 
including the on-site nonpolynomial nonlinearity and 
long-range DD interactions is formulated in Section II. 
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Results for on-site and inter-site solitons, including the 
study of their stability, and the possibility of the col- 
lapse onset, are presented in Section III. The mobility of 
discrete solitons in the DNPSE model is also briefly con- 
sidered in Section III. The paper is concluded by Section 
IV. 



II. THE MODEL 

Adding the long-range DD interaction to the scaled 
nonpolynomial Schrodinger equation (in its continual 
form Jy, 01) leads to the following ID equation: 



the repulsive DD interactions). Two dynamical quanti- 
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Here, Vq and t: / q are the strength and period of the lon- 
gitudinal optical lattice potential, F{z,t) = ^/\j\f{z,t), 
where / is the ID mean-field wave function subject 
to normalization \f{z,t)\'^ dz = 1, and 7 — 



—2Nas \/ muj± / h is the effective strength of the local 
interaction, with N the total number of atoms in the 
condensate, the scattering length of atomic collisions 
{ttg < corresponds to the attraction), m the atom mass, 
and Lj± the transverse trapping frequency @. Further, 
H = sgn(as) is the sign of the local interaction, and 
G = g (1 — 3cos^ 6) is the coefficient which defines the 
ratio of the DD and contact interactions, where 17 is a 
positive coefficient and 6 the angle between the z axis 
and the orientation of the dipoles. One obvious case of 
interest in the ID geometry corresponds to the dipoles 
polarized along the z axis, i.e., 6 ~ and G = —2g, when 
the long-range interaction is attractive. Another inter- 
esting case corresponds to the orientation of the dipoles 
perpendicular to the z axis (0 = 7r/2 , i.e. G ~ g), which 
implies the repulsive DD interaction. 

Assuming a sufficiently deep optical lattice potential, 
and approximating the wave function by a superposition 
of localized Wannicr modes, similar to how it was done 
in Refs. J)] and one can derive the discrete version 
of Eq. H]), i.e., the DNPSE with the DD term: 
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where the linear-coupling G and the ratio between DD 
and contact interaction coefficients F — Gj I7I, are de- 
fined as in Ref. [21J (F > for attractive and F < for 



ties are conserved by Eq. ([2]), ot2., 
and Hamiltonian 
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Note that the staggering transformation, Fn = 
(— 1)" exp (— 4iCi) can be used to change the sign 
of C if it is negative, but this transformation cannot be 
used to invert the signs of nonlinearity coefficients in Eq. 
0- 

Experimentally adjustable parameters are the rela- 
tive strength of the DD/contact interactions, F, and the 
norm of stationary wave function, P. The latter may 
be expressed in terms of the total number of atoms in 
the condensate, [i, [2l|: P ~ N \as\ ^/muiT/h. In 
particular, « 5 nm for ^^Cr atoms (far from the 
FR); assuming the transverse-confinement width to be 
{y^muj±/h)^^ 5 /im, we conclude that P ~ 1 may cor- 
respond to ^ 1000 atoms in the condensate. A typical 
value of the relative interaction strength, if estimated as 
the ratio of the effective scattering lengths corresponding 
to the DD and contact interactions in the ^^Cr conden- 
sate, without the use of the FR technique, is |F| ~ 0.15 
[lll |. Actually, F can be made both positive and negative, 
and its absolute value may be altered within broad limits 
by means of the FR [l3|. In particular, |F| can be given 
very large values in the experimentally possible situation 
when the strength of the contact interactions is almost 
nuUified with the help of the FR 

Stationary solutions to Eq. ([J), with chemical poten- 
tial /i, are sought for as P„ = u„ exp(— i/if), with real 
discrete function u„ satisfying a stationary equation. 
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In the case of the attractive contact interaction (K = +1), 
ui cannot exceed the maximum value, (u?,) = 1, as 

n ' V / max ' 

seen from Eq. ([4]). In fact, the presence of the singularity 
in Eq. (H]) at |P„| =1 makes it possible to study the 
onset of collapse in the framework of this equation fo'] . 

Stationary equation ([4]) was numerically solved by 
an algorithm based on the modified Powell minimiza- 
tion method 0, [2ll |. Initial ansdtze used to con- 
struct on-site and inter-site-centered discrete solitons 
were taken, respectively, as = (..., 0, A, 0, ...) 

and (..., 0, A, A, 0, ...), where A is a real constant ob- 
tained from Eq. ^ in the corresponding approximation. 
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Results reported below were obtained in the lattice com- 
posed of 101 or 100 sites, for the on-site and inter-site 
configurations, respectively. It was checked that the re- 
sults do not alter if a larger lattice had been used. 

III. RESULTS AND DISCUSSION 
A. The case of the attractive contact interaction 

Families of fundamental unstaggered solitons of on-site 
and inter-site types for the local attraction (K = +1) 
and cither sign of the DD interaction were obtained from 
the numerical solution of Eq. In contrast to the 

ID discrete model with the cubic on-site nonlinearity, 
where bright unstaggered solitons exist only for suffi- 
ciently weak repulsive DD interaction (2ll | , in the present 
model solitons can be also found if the repulsive DD inter- 
action is strong. The solitons were categorized as stable 
ones if they met two conditions: the slope ( Vakhitov- 
Kolokolov) criterion, according to which the slope of the 
P(/i) dependence must be negative (or may be very close 
to zero, see below), dP/d/j, < 0, and, simultaneously, the 
spectral condition, according to which the correspond- 
ing eigenvalues, found from linearized equations for small 
perturbations, must not have a positive real part (2ll.l23l|. 

In order to draw general conclusions about the influ- 
ence of the DD interaction on the solitons, the analysis 
is presented below for two different values of the lattice 
couphng constant C, viz., C = 0.8 and C = 0.2. These 
cases correspond, respectively, to the proximity to the 
continuum limit, and to a strongly discrete system. 

1. The quasi- continual model (C = 0.8^ 

A global characteristic of soliton families is the /'(/i) 
dependence, i.e., the scaled norm as a function of the 
chemical potential. For C = 0.8 and four different val- 
ues of DD parameter F, the P{ii) curves are displayed in 
Fig. [TJ In the absence of the DD interaction, F = 0, two 
subfamilies of on-site solitons are found, one (which oc- 
cupies a narrow interval of /x) obeying the slope criterion, 
and the other one violating it, in a broad interval. With 
the increase of the strength of the attractive DD interac- 
tions (F = 5), the region where the slope condition is met 
spreads out, and, when the DD interaction is dominant 
(F = 12), the condition is satisfied for all on-site soli- 
tons. On the other hand, the slope condition is fulfilled 
for all inter-site solitons, regardless of the value of F. It 
is worthy to note that the difference between the P(/i) 
curves for the on-site and inter-site solitons vanishes as 
the attractive DD interaction strengthens. 

In the case of the repulsive DD interaction, the P{n) 
curves for on-site and inter-site solitons are completely 
separated, see Fig. [Ud). Actually, the slope of P{^j) 
curves for both on-site and inter-site solitons tends to be 
very small in this case, which makes the application of the 



slope criterion doubtful. In these areas, the solitons fea- 
ture the amplitude close to the upper limit admitted by 
the model, w^ax = 1? ^'^'^ ^ very small width, correspond- 
ing to a situation when nearly all atoms are collected in 
a single well of the underlying potential lattice. 

The spectral stability was examined by linearizing Eq. 
^ for small perturbations 5Fn around the soliton, and 
finding the respective eigenvalues in a numerical form. 
The results of the stability analysis for the on-site soli- 
tons, with C = 0.8, are summarized in Fig. [2l in the 
form of the stability diagram in the plane of (/z, F), where 
contours of constant norm P for the on-site solitons are 
included too. The collapse condition, u^ax = 1) is at- 
tained at the black solid line, which is, simultaneously, 
a stability border. In direct simulations, the on-site soli- 
tons which are predicted to be stable survived as long as 
the simulations ran [Fig. El^a)], while the solitons clas- 
sified as unstable ones suffered the collapse (destruction 
of the solution after attaining the level of u^ax = 1) in a 
finite time. 

On the other hand, all inter-site solitons are unstable, 
as the spectrum of eigenvalues for small perturbations 
around them always contains real eigenvalues. How- 
ever, unstable inter-site solitons which have stable on- 
site counterparts with the same norm avoid the collapse, 
evolving into robust breathers oscillating around the cor- 
responding stable on-site solitons, as shown in Fig. [31(b). 
Unstable inter-site solitons do collapse if no stable on- 
site soliton with the same norm can be found, see Fig. 
[3fc). Therefore, the border line for the collapse of the 
inter-site solitons coincides with the stability border for 
on-site solitons in Fig. [51 

For the repulsive DD interaction, the spectral stabil- 
ity condition for all inter-site solitons does not hold ei- 
ther. Because, in this case, curves P(^) for the inter-site 
and on-site solitons are strongly separated [Fig[TKd)], un- 
stable inter-site solitons do not find stable on-site coun- 
terparts with the same norm, therefore they suffer the 
collapse. As for the on-site solitons, there exists a re- 
gion where the spectral stability condition holds for them. 
This region expands as the repulsive DD interaction gets 
stronger, although the respective curve P{^) becomes al- 
most horizontal, and the slope condition cannot be ac- 
curately verified. Direct numerical simulations confirm 
the predictions of the stability analysis for the on-site 
solitons in this case too. 



2. The strongly discrete model (C — 0.2) 

In this case, the slope condition is satisfied for all on- 
site solitons in the absence of the DD interaction (F = 0), 
see Fig. [H^a). As the strength of the attractive DD in- 
teraction grows, a pair of subfamilies emerge, the slope- 
stable and unstable ones, the respective P{n) curves be- 
ing similar to those observed in the case of the strong 
couphng, C = 0.8 [Fig. [HJb)]. Eventually, when the 
attractive DD interaction becomes dominant, the region 
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where the slope condition is satisfied spreads over the 
entire parameter space, as seen in Fig. Hfc). 

In the model with C — 0.2, all inter-site solitons again 
satisfy the slope condition. With the strengthening of the 
DD interaction, the separation between the P(/i) curves 
for the on-site branch, which satisfies the slope condi- 
tion, and its inter-site counterpart vanishes. On the other 
hand, in the case of the repulsive DD interaction, the cor- 
responding P(/i) curves for on-site and inter-site soliton 
families are strongly separated, see Fig. [DJd). 

Results of the stability analysis for the on-site soli- 
tons are summarized in the stability diagram displayed 
in Fig. Efa) in parametric space (r,/i). Unlike the case 
of C = 0.8, cf. Fig. [21 in the present case there ap- 
pears a region where the spectral stability condition is 
violated for on-site solitons as the DD interaction grows 
stronger, as well as a region where the stability condi- 
tion holds for inter-site solitons. These regions disappear 
again with the further growth of F. Also in contrast to 
the case of C = 0.8, the line at which the collapse is at- 
tained in the family of on-site solitons does not coincide 
with the border between stable and unstable parts of the 
family. Rather, the collapse line passes trough the un- 
stable region where the spectral-stability condition does 
not hold, see Fig. [Sja). Direct simulations demonstrate 
that unstable on-site solitons which have stable inter-site 
counterparts with the equal norm evolve into breathers 
oscillating around the stable counterparts. On the other 
hand, if unstable on-site solitons cannot find stable inter- 
site counterparts with the same (or close) value of the 
norm, they undergo the collapse. In fact, the existence of 
the two different scenarios of the instability development 
- the formation of the breather and collapse - explains 
the above-mentioned fact that the collapse line does not 
coincide with the instability border. In the present case 
too, direct simulations corroborate the stability of those 
solitons which do not have unstable eigenvalues. 

Extending the above-mentioned trend, unstable inter- 
site solitons which have stable on-site counterparts with 
the same norm evolve into persistent breathers, while 
those unstable solitons that are devoid of stable equal- 
norm counterparts suffer the collapse. Accordingly, the 
line at which the collapse is attained does not coincide 
with the instability border, as seen in Fig. [DJb). It is 
worthy to note the existence of a stability region for inter- 
site solitons in Fig. [DJb) which is adjacent to the collapse 
line. In the latter case, the stable inter-site solitons do 
not have on-site counterparts with the same value of the 
norm. 

For the repulsive DD interactions, the spectral stability 
condition always holds for on-site solitons and does not 
hold for inter-site modes, similar to the case of C = 0.8. 
Again, in some part of the parameter space, the corre- 
sponding P{n) curves for the on-site solitons are nearly 
horizontal lines with zero slope, being completely sepa- 
rated from the -P(/x) line for the inter-site modes. Direct 
simulations demonstrate that the on-site solitons are in- 
deed stable in this case, while the unstable inter-site soli- 



tons undergo the collapse. 



B. The case of repulsive contact interactions 

If the local interaction is repulsive, the existence of 
the unstaggered solitons may only be supported by the 
attractive DD interaction. In the model with the self- 
repulsive cubic on-site nonlinearity, unstaggered solitons 
were found only when the relative strength of the attrac- 
tive DD interaction was large enough, F > 0.4 [2TI|. In the 
present model, unstaggered solitons (with large ampli- 
tudes) are obtained for smaller values of F as well. Nev- 
ertheless, general results obtained in the present model 
for the case of the local repulsion are not drastically dif- 
ferent from those reported in the model with the cubic 
on-site nonlinearity in Ref. (2l| . In particular, differ- 
ences between P(/Lt) curves for on-site and inter-site soli- 
tons are small at all values of /i, the slope condition is 
always satisfied for both types of the solitons, and there 
is an exchange between regions where the spectral stabil- 
ity condition is fulfilled for on-site and inter-site solitons. 
Further, unstable solitons evolve into breathers oscillat- 
ing around their stable counterparts. The similarity of 
these results for discrete solitons in the present DNPSE 
model and its cubic counterpart is not surprising, as in 
the case of the local repulsion (unlike attraction) there is 
no dramatic difference between the nonpolynomial and 
cubic nonlinearities, therefore the competition of the lo- 
cal term of either type (nonpolynomial or cubic) with the 
DD attraction gives rise to similar solitary modes. 



C. Moving discrete solitons 

Finally, we briefly summarize results concerning mo- 
bility of localized modes, with respect to the concept of 
the Peierls-Nabarro barrier [23| . To this end, we follow 
the lines of the analysis developed in Refs. Q and 
Examination of the mobility has shown that the Peierls- 
Nabarro barrier in the DNPSE model depends on the 
strength of the DD interaction in the same way as it was 
in the discrete model with the cubic on-site nonlinearity. 
Namely, in the case of the local attraction, the repul- 
sive or attractive DD interaction decreases or increases, 
respectively, a region in plane (P, /i) where mobile local- 
ized modes can be found. On the other hand, in the case 
of the contact repulsion, all localized modes can be set in 
motion by an initial kick, regardless of the value of DD 
coefficient F. In all cases when mobile discrete solitons 
exist, the vanishing of the Peierls-Nabarro barrier coin- 
cides with the disappearance of the separation between 
P(/z) curves for the on-site and inter-site soliton families, 

cf. Refs. d,!!!!. 
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IV. CONCLUSION 

The purpose of this work was to achieve a better 

understanding on the influence of the long-range DD 
(dipole-dipole) interactions on the stability and collapse 
of localized nonlinear modes in the cigar-shaped Bose- 
Einstein condensate trapped in a deep optical-lattice po- 
tential. To this end, wo have introduced the model based 
on the one-dimonsional DNPSE (discrete nonpolynomial 
Schrodinger equation), which includes the contact (on- 



site) and DD nonlinear terms. Both attractive and re- 
pulsive signs of the contact and DD interactions were 
considered. The main conclusion is that the presence of 
the attractive DD interaction enhances the soliton's sta- 
bility and helps to prevent the collapse. Our analysis 
was limited to unstaggered solitons, the consideration of 
staggered localized modes being a subject of a separate 
work. 

G.G., A.M. and Lj.H. acknowledge support from the 
Ministry of Science, Serbia (through project 141034). 
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Figures 




FIG. 1: The P{fi) (norm versus chemical potential) depen- 
dencies for families of on- and inter-site unstaggered solitons 
(the solid and dashed lines, respectively) in the case of the at- 
tractive contact interaction, for C = 0.8 and relative strength 
of the DD interaction T = (a), 5 (b), 12 (c), and -5 (d). 
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FIG. 2: The stability /collapse diagram for on-site solitons, as 
determined by the stability analysis performed for C — 0.8. 
The gray and white areas are occupied by unstable and stable 
solitons, respectively. The soliton norm keeps constant values 
along thin contour lines, as indicated in the figure. The chain 
of bold dots connected by the black line shows a curve along 
which the collapse is attained. 




FIG. 3: (a) An example of a stable on-site soliton. (b,c) 
Development of the instability of inter-site solitons: (b) the 
case when a stable on-site soliton exists whose norm is equal 
to that of the unstable inter-site sooliton. In this situation, 
the unstable soliton evolves into a breather oscillating around 
the stable on-site soliton. (c) The case without the stable on- 
site counterpart with the equal norm. In the latter case, the 
unstable soliton collapses. 
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FIG. 4: P(/i) dependencies for families of on- and inter-site 
(solid and dashed lines, respectively) unstaggered solitons in 
the strongly discrete model (C — 0.2) with the attractive con- 
tact interaction. The relative strength of the DD interaction 
is r = (a), 5 (b), 15 (c), and -5 (d). 
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FIG. 5: The stability/collapse diagram for (a) on-site (a) and 
inter-site (b) unstaggered solitons in the strongly discrete ver- 
sion of the model, with C = 0.2. The notation has the same 
meaning as in Fig. [5] 



